Semiclassical methods for multi-dimensional systems bounded by finite potentials 
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This work studies the semiclassical methods in multi-dimensional quantum systems bounded 
by finite potentials. By replacing the Maslov index by the scattering phase, the modified transfer 
operator method gives rather accurate corrections to the quantum energies of the circular and square 
potential pots of finite heights. The result justifies the proposed scattering phase correction which 
paves the way for correcting other semiclassical methods based on Green functions, like Gutzwiller 
trace formula, dynamical zeta functions, and Landauer-Biittiker formula. 
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Semiclassical approaches are techniques to study quan- 
tum systems in classical limit Q. In these approaches, 
physical quantities are expressed by the Green function 
of the system and replaced by its underlying classical 
trajectories under stationary phase approximation. Well 
known examples range over Gutzwiller trace formula, dy- 
namical zeta function, and transfer operators, which are 
aimed at determining the quantized energies of closed 
systems Q. For open systems, the most prominent exam- 
ple might be the Landauer-Biitticker formula for charge 
current transport through open quantum dots 0. Re- 
cently semiclassical approaches also have been used to 
study spin current transport and spin dynamics, which 
clarifies the suppression of D'yakonov-PereF spin relax- 
ation in mesoscopic quantum dots Q. All these ap- 
proaches can relate quantum problems to the ergodicity 
property of their corresponding classical dynamics and 
gives more transparent pictures to complexphenomena, 
including the signatures of quantum chaos [ij . 

Besides this conceptual contribution, some semiclas- 
sical methods provide efhcient numerical techniques for 
less time-consuming calculations. Their results are es- 
pecially accurate in mesoscopic systems in which the de 
Broglie wavelength of the Fermi electrons is much shorter 
the sample size. The devices in this scale are often fab- 
ricated by lithography or controlled by confining poten- 
tials. A convenient theoretical approach to study these 
quantum systems is assuming them to be bounded by 
infinite potential. This largely simplifies the formulism 
of the semiclassical methods, since the quantum particle 
has the same phase change for arbitrary particle energy 
after it is reflected by this potential. This phase change is 
carried in the well known Maslov index. However, for real 
systems beyond this assumption, the accuracy of the con- 
ventional semiclassical methods using Maslov index be- 
come out of control. The usual way is going back to solve 
the Schrodinger equation with a boundary combined with 
Dirichlct and von Ncuman conditions. But one could ask 
whether the conventional semiclassical methods still work 
after taking certain effective correction. 



A natural candidate for this correction is an effective 
phase change in the wave function after it is bounced 
back by the potential. This phase can be understood as 
the dwell time of the quantum particle penetrating into 
and staying inside the potential barrier. Indeed, in the 
one-dimensional (ID) potential well, it has been shown 
that replacing the Maslov index in the Green function by 
a scattering phase, this function gives an exact quanti- 
zation rule identical with the Wentzel-Kramers-Brillouin 
(WKB) method 0. This relation stirs up the motivation 
how to extend this result to multi-dimensional systems. 
The current paper demonstrates this extension in the ex- 
ample of Bogomolny's transfer operator (BTO) method 
and tested it in a circular billiard and a square billiard 
bounded by potential pots of finite hight. The calculated 
energies are surprisingly accurate, which justifies the sug- 
gested scattering phase correction for multi-dimensional 
systems. Although the result is presented in the BTO 
method, it should be valid for all semiclassical methods 
as long as they are based on the Green function. 

The WKB method might be the simplest semiclassical 
method. For a particle of mass m and energy E bounded 
by an ID potential V{x), the solution of its Schrodinger 
equation can be approximated by the WBK wave func- 
tion [3, 
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with p{x) ~ \/2m\E ■ 
length A(x) 
potential. 



y(a;)], if the de Broglie wave- 
2T:h/p{x) varies slowly compared to the 
This requirement is usually violated at the 
classical turning point a;o i where the momentum changes 
sign. At that point one has E — V{xq) which gives rise 
to a vanishing p{xq) and a singular function 'il'{x). If the 
potential varies slowly around xq, the exponentially de- 
creasing real wave function outside this point should be 
associated with the oscillating wave function inside this 
point. This enforces the incident wave function to take 
a scattering phase after reflection. This phase is equal 
to 7r/2 in the semiclassical (short wave) limit If the 
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particle is reflected back and forth between two turning 
points xi and X2 on two potential barriers, the particle 
should take two scattering phases 0i = tt/2 during 
one period of motion. The total phase then equals 2mr 
with an integer n, which leads to the WKB quantization 
condition 
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where the Maslov index — 2 corresponds to the two 
reflections during one period of particle motion. 

If the potential does not vary sufficiently slowly at the 
turning points, for example the step function barrier, the 
phase change 7r/2 is no longer a good approximation. A 
general scattering phase should be determined quantum- 
mechanically [^, 
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where the phase changes </>! {E) and 02 (E) at the turning 
points xi respectively X2 become a function of E. As an 
example, let this particle move in an ID finite square well 
with V{x) = for < x < L and Vq otherwise, where L is 
the well width and Vq > is the potential height. Solving 
the Schrodinger equation with the continuous boundary 
condition, the scattering phase can be calculated. 
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Substituting Eq. (gj) into (pi{E) and (p2{E) in Eq. 
one obtains an exact quantization rule for the ID finite 
square well. 

For general fc-dimensional systems with k > 2, many 
quantization rules have been derived to extend the WKB 
method, including Gutzwiller trace formula, dynamical 
zeta function, and transfer operators All of them are 
based on the semiclassical Green function 



G{r,r';E) = 
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which is a sum over all classical trajectories 7 starting 
from r' and ending at r 0,0. The function Sc\{y,y' ; E) 
therein is the action of the particle from r' to r along 
7. The Maslov index /i stands for the total number of 
the turning points between r' and r along 7. For ID 
finite square well, the quantization rule derived from Eq. 
(O with the replacement /x^ by in Eq. igj is identical 
with the quantization rule of Eq. Q with the same phase 
correction (j)s 

For multi-dimensional systems, intuitively, the magni- 
tude of the scattering phase should depend on the inci- 
dent angle of the particle wave. Naively one could suggest 



the phase 0s in Eq. @ to be only related to the perpen- 
dicular component of the incident wave with respect to 
the boundary. Figure ^a) shows an example of a 2D 
step potential Vc{x,y) = VqQ{x) with the potential high 
Vo > for a; > and for cc < 0, where Q{x) is the 
Heviside function. The momentum px perpendicular to 
the boundary has energy Ep — p'^/2m, which will replace 
the total energy E in Eq. for multi-dimensional sys- 
tems. This phase is a value between and tt, depending 
on the ratio Ep/Vo, as shown in Fig. ^b). If the po- 
tential barrier is high, that is Ep/Vo Ij the scattering 
phase will approach tt, which is the same as the infinitely 
high potential barrier. 
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FIG. 1; (a) The incident state \px,Py) and the reflected state 
I — PxjPv) e"'^" on the confining potential Vo- (b) The scat- 
tering phase (j)s of the reflected state is a function of Ep/Vo- 



Next, the scattering phase for the perpendicular wave 
component discussed above will be merged into BTO 
method. Let us consider again the particle of mass m 
and energy E moving in a fc-dimensional system. Se- 
lect a Poincare section (PS) S in the configuration space 
of this system, such that almost all classical trajecto- 
ries pass this section lOj. The original transfer operator 
T{E) is defined as the integral operator 



T{E)^P{q)^ / T{q,q';E)i;{q')dq', 
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acting on some function 'ip{q') on S. The integral kernel, 
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is defined as the sum over all possible classical trajectories 
7's from the initial point g' 6 S to the final point g 6 S. 
The action Sci{q,q';E) is the same as before in Eq. 
The Maslov index v counts the number of the crossing 
points of 7 through S from the same side of S. According 
to the BTO method, the zeros of the Fredholm determi- 
nant I det(l — T{E))\ of the transfer operator T{E) are 
the energy eigenvalues of the quantum system. 
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The kernel in Eq. ^ is derived from the semiclas- 
sical Green function in Eq. (0), which can be divided 
into two parts, G{q,q'\E) = G°^''iq, q'; E) -Go{q,q';E), 
where G°^'^{q, q'; E) is the contribution from long trajec- 
tories and Go((7, q'] E) is the contribution form short tra- 
jectories. After coordinate deduction from /c-dimensional 
space to (fc — l)-dimensional space on E 0, the quanti- 
zation condition | dei{G{q,q' ] E))\ = can be expressed 
as I det(l — T{E))\ — 0, where the identity operator 
comes from Gq and the operator T{E) originates from 
Qosc_ The entire derivation holds the same when the 
phase in G{q,q';E) is modified. The transfer operator 
Tm (E) modified with the scattering phase for perpendic- 
ular wave component then has the kernel 
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The efficiency of this modified operator will be tested in 
the following two 2D integrable systems. 
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FIG. 2: The upper part shows the Fredholm determinant 
I det(l — Tm{E))\ of the modified transfer operator Tm{E) for 
the potential pot in the inset for Vb — 20, 50, and oo with 
R — 2. The points on the three bottom lines depict the exact 
quantum energies for Vb = 20, 50, and oo with R — 2. 



The first system is a 2D circular quantum dot bounded 
by a finite potential, as shown in the inset of Fig. |21 The 
potential pot has radius R and height Vb, that is V{r) = 
Vb for r > R and V{r) = for r < R. Analytically 
the energy eigenvalues of its Schrodinger equation can 
be calculated by matching the boundary conditions at 
radius r — R Setting Planck constant h = 1, the 

mass m — 1, and the pot radius i? = 2, the exact energy 
eigenvalues for different potential height Vb — 20, 50, 
and oo are denoted on the three bottom fines of Fig. |21 
The dotted, dashed, and solid curves in the upper part 



of Fig. 12 are the Fredholm determinant | det(l — %n{E))] 
of the modified transfer operator for Vb = 20, 50, and oo. 
The zeros of these functions determine the semiclassical 
quantum energies of the systems. 
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FIG. 3: The upper part shows the Fredholm determinant 
det(l — Tm{E))\ of the modified transfer operator Tm{E) for 
the square potential pot in the inset with the potential con- 
figuration {V^,Vy) = (20,20), (50,50), (20, oo), and (oo,oo) 
where Vxy = Vx + Vy. The points on the four bottom lines 
depict the exact quantum energies for {Vx,Vy) = (20,20), 
(50,50), (20, c»), and (cx),cx)). 



The second system is the 2D square quantum dot with 
the confining potential as shown in the inset of Fig. O 
Therein the xy-plane is separated into nine regions with 
different constant potential heights Vx, Vy, and Vxy at 
each region, where Vx > 0, Vy > 0, and Vxy = Vx + Vy. 
This 2D problem can be reduced to two independent 
ID finite wells with potential heights Vx and Vy and 
solved separately Combining the eigenvalues of 

these two separated systems, the total quantum ener- 
gies of this 2D system for iVx,Vy) = (20,20), (50,50), 
(20, oo), and (oo, oo) are determined and depicted on the 
four bottom lines in Fig. |3| where 7i = m = 1 as be- 
fore and the well length L is normalized by the condition 
2tt^ /L"^ ~ 1. The upper part of this Figure shows the 
Fredholm determinants | det(l — %n{E))\ of the modified 
BTO. The dashed, dash-dotted, dotted, and solid curves 
represent these functions for the potential configurations 
{Vx,Vy) = (20,20), (50,50), (oo,oo), and (20, oo). 

Figure 2 and 3 shows surprisingly accurate quantum 
energies after the scattering phase correction. Taking 
Fig. 2 as example, the exact 10-th energy has a remark- 
able left shift after the potential height is reduced from 
Vb = oo to Vb = 20. This shows how large the error could 
be when simplifying a deep potential pot, even with the 
large ratio Vq/R— 10, by an infinitely high pot. Without 
scattering phase correction, the semiclassically calculated 
energies are the zeros of the solid curve. The 10-th zero 
of this curve is close to the 10-th point on the bottom 
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line of Vo = oo, but far apart from the 10-th point for 
Vq = 20. However, after taking the scattering phase, the 
determinant (dotted curve) shifts leftward quite a lot and 
its zeros largely approach the exact energies for Vq = 20. 

Quantitatively we can define values to characterize 
these errors. Suppose Ei is the i-th exact energy of a 
quantum system bounded by some potential pot of height 
Vb < oo and Ei is the corresponding semiclassical en- 
ergy approximated by the modified BTO. Let the value 
Si = Ei~ Ei be the difference between these two energies 
and the ratio = Si/Ei is the error of the i-th eigen- 
value. For the special case of Vb = oo the values defined 
above are furnished with a superscript oo as E°°, E°°, 
and 6°° = E°° — E°° respectively. The relative error Ti 
of the i-th energy is then defined as the ratio 



\E^ - E? 



(9) 



The denominator denotes the exact energy shift after the 
potential height is reduced from oo to Vq. The numerator 
represents the error 6i for Vb < oo subtracted by the basic 
semiclassical error S°° from infinite potential. 
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FIG. 4: The dotted and dashed curves in the main plot de- 
note the relative errors Ti of the 2D circular potential pots 
with potential height 50 and 20, as shown in the inset of Fig. 
13 The dash-dotted and solid curves represent Fi for the 2D 
square potential pots of configuration {Vx, Vy) = (50, 50) and 
(20,20), as shown in the inset of Fig. |21 The corresponding 
errors Si of these systems are shown in the inset. 



According to these definitions, the errors 6i of the low- 
est ten energies from Fig. |21 and 13 are plotted in the 
inset of Fig. ^ and their corresponding relative errors F^ 
are plotted in its main figure. All 6i are bounded by 5% 
and all F^ are bounded by 20%. Roughly speaking, the 
scattering phase in the modified BTO has corrected at 
least 80% of the energy error due to the potential reduc- 
tion from oo to Vb. This justifies the scattering phase 
correction proposed above for finitely confined systems. 



Notably, one cannot expect 100% correction in this 
first order correction. Remember that there still exists 
other degrees of freedom in the system, which are not 
included in the scattering phase. For instance, the mod- 
ified transfer operator T„i{E) is the same for square po- 
tential pots of different Vxy, although they have different 
quantum energies. This Vxy difference can only be distin- 
guished from higher order corrections beyond the current 
scattering phase. Furthermore, if the particle energy E 
is close to the potential height Vb, the particle can pen- 
etrate into the potential well quite long and the wave 
property of the particle prevails its particle property. In 
this regime, the correction deviation from the modified 
BTO increases. However, that is not because this first 
correction is wrong, but because higher order corrections 
are required in the semiclassical approach. 

Finally, almost all semiclassical methods are based on 
the Green function of the quantum system. The success- 
ful result in the modified BOT gives a clear direction for 
extending the quantum correction to other semiclassical 
methods, including Gutzwiller trace formula, dynamical 
zeta functions, and Landauer-Biittiker formula. 
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